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Steepening mass profiles, dark matter and environment of 
X-ray bright elliptical galaxies 
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ABSTRACT 

We use a new non-parametric Bayesian approach to obtain the most probable mass 
distributions and circular velocity curves along with their confidence ranges, given 
deprojected density and temperature profiles of the hot gas surrounding X-ray bright 
elliptical galaxies. For a sample of six X-ray bright ellipticals, we find that all cir- 
cular velocity curves are rising in the outer parts due to a combination of a rising 
temperature profile and a logarithmic pressure gradient that increases in magnitude. 
Therefore at large radii, mass density profiles rise more steeply than isothermal pro- 
files, implying that we are probing the more massive group-sized haloes in which these 
galaxies are embedded. Comparing the circular velocity curves we obtain from X-rays 
to those obtained from dynamical models, we find that the former are often lower 
in the central ~ 10 kpc. This is probably due to a combination of: i) Non-thermal 
contributions of up to ~ 35% in the pressure (with stronger effects in NGC 4486), 
ii) multiple-temperature components in the hot gas, iii) incomplete kinematic spatial 
coverage in the dynamical models, and iv) mass profiles that are insufficiently general 
in the dynamical modelling. Complementing the total mass information from the X- 
rays with photometry and stellar population models to infer the dark matter content, 
we find evidence for massive dark matter haloes with dark matter mass fractions of 
~ 35-80% at 2R e , rising to a maximum of 80-90% at the outermost radii. We also 
find that the six galaxies follow a Tully-Fisher relation with slope ~ 4 and that their 
circular velocities at lR e correlate strongly with the velocity dispersion of the local 
environment. As a result, the galaxy luminosity at lR e also correlates with the ve- 
locity dispersion of the environment. These relations suggest a close link between the 
properties of central X-ray bright elliptical galaxies and their environments. 

Key words: galaxies: elliptical and lenticular, CD - galaxies: evolution - X-rays: 
galaxies - galaxies: stellar content - dark matter 



1 INTRODUCTION 

X-ray bright elliptical galaxies are massive galaxies, thought 
to be among the most evolved systems in our Universe with 
a complex formation history. They are believed to grow from 
mergers between smaller galaxies, and after star formation 
ceases at z ~ 1-2, they become larger and less compact 
through the significant accretion of stellar material from 
neighbouring sm aller systems (e.g. Ivan Dokkum et al.ll200a ; 
iNaab et al.ll2009l ). 

Knowledge of the total mass distributions of galaxies 
gives us insight into their formation history at several dif- 
ferent levels. At the most basic level, it tells us the combined 
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mass of these collapsed systems. If we incorporate informa- 
tion from photometry and stellar population models, we can 
disentangle the respective luminous and dark matter com- 
ponents. Comparing their properties with what is expected 
from simulations will help place constraints on current theo- 
ries of galaxy evolution. The mass profile could also be used 
as input in dynamical models of galaxies, therefore miti- 
gating the usual mass-anisotropy degeneracy. This would 
provide more stringent constraints on the orbital structure, 
which serves as an imprint of the processes that occurred 
in the past to create the galaxy. Finally, it will give us an 
insight into the relations between the masses of collapsed 
systems and the environments they reside in. 

There are several methods in the literature for ob- 
taining the mass distributions of elliptical galaxies. Dy- 
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Galaxy 
(1) 



Classification 

(2) 



Distance (Mpc) 
(3) 



logLx 

(4) 



ffcnv Source 
(5) (6) 



NGC 1399 


cD; El pec 


NGC 1407 


E0 


NGC 4472 


E2 


NGC 4486 


E0-E1 pec 


NGC 4649 


E2 


NGC 5846 


E0-E1 



19.95 
28.84 
16.29 
16.07 
16.83 
24.20 



41.63 


370 


D01 


41.00 


387 


T06 


41.43 


607 


G99 


42.95 


762 


G99 


41.28 


702 


G99 


41.65 


322 


M05 



Tabl e 1. S a mple of X-ray bright elliptical galaxie s: (1) Galaxy name, (2 ) galaxy classification from NED, (3) distance taken from 
iTonrv et al] l|200lh . (4) X-ray luminosity taken from lO'Sullivan et al . (2001), (5) veloc ity dis persion of the surround ing en vironment and 

(6) so urce of the velocity dispersio n values. References are D01 JDrinkwater et al,ll200lT ). T06 JTrentham et alj2006l ). G99 I IGavazzi et al] 
1 19991 ). M05 JMahdavi et alj|2005t) . 
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namical models can be constructed by su perposing a li- 
brary of orbit s (e.g. iRix et al.lll997l; iGebhardt et alj|2003t 



I200J. 
g. iDeionghe et all 1 19961 ; iGerhard et al] iFukazawa et al.l | 
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hydrostatic equilibrium using a procedure that is as free 
fro m systematic bias es as possible. 



or distribu- 



1 19981 ; iKronawitter et al.l 20001). or by construc ting a sys 



Cowie et al] (1 19871) 

' 20061 ); iNagino fc Matsushital 
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tern of particles (|de Lorenzi et al] I200S1 . 120091 ) such that 
the projection of the system best reproduces observed 
surface-brightness and kinematic profiles. These models 
give the mass distribution and orbital structure of the 
galaxies simultaneously. Strong lensing gives the pro- 
jected mass within the Einstein ring and weak lensing 
studies provi de mass density profile s for a sample of 
galax i es (e.g. | Treu fc Koopmansl |2004|; iMandelbaum et al] 



I2006J ; iGavazzi et al] 120071 ; iKoopmans et al] I2009T I. The 
properties of X-ray bright elliptical galaxies allow an 
additional possibility for obtai ning mass distributions 



(e.g. INulsen fc Bohringerl Il995l; IFukazawa et al] 120061 : 



IChurazov et al.l l200Sl : INagino fc Matsushital 120091 ). The X- 
ray spectra are dominated by lines and by continuous emis- 
sion from thermal bremsstrahlung radiation produced in the 
surrounding halo of hot gas. The spectra can be deprojected 
and fitted to derive 3-D temperature and density profiles of 
the gas, and if we assume it is in hydrostatic equilibrium, 
we can derive the enclosed mass profile. 

Generally the literature points towards isothermal mass 
distributions in elliptical galaxies, with a conspiracy be- 
tween the luminous and dark matter components resulting 
in flat circular v elocity curves similar to t hat found in spi- 
ral g a laxies (e.g . Krona witter et a l. 2000; Koopman s et al] 
120091 ; IChurazov et"afl l201Ch . Fukazaw aetTafl l|2006D how- 



ever, determined mass distributions of a sample of 53 el- 
liptical galaxies from Chandra and XMM-Newton observa- 
tions, and found that they are better described by power 
laws with an index > 1, i.e. their circular velocity curves 
range bet ween flat and rising. The part icle-based dynamical 
models of lde Lorenzi et al] ([2008, 2009) found falling circu- 
lar velocity curves for the in termediate- lum i nosity galaxies, 
NGC 3379 and NGC 4697. IThomas et all (|2007l ) found a 
range of slopes in the elliptical galaxies belonging to the 
Coma cluster, from their orbit-based dynamical models. 

In this paper, we examine the isothermality of mass 
distributions from X-ray observations and whether they are 
accurate enough to derive dark matter mass profiles, use 
in dynamical modelling, and determine global properties of 
elliptical galaxies. To address these issues we need to apply 



used methods where the measured temperature and density 
profiles are parameterised before applying hydrostatic 
equilibrium, to circumvent differ entiating the o b served 
profiles, which are often noisy. IChurazov et al] 12008) 
avoided differentia tion by inte rpolating to obt ain only the 
potential profile. INulsen fc Bohringerl |l995j) found the 
most likely mass profile of NGC 4486 (M87) from temper- 
atur e and density pro f iles us ing a non-parametric method, 
and Irlumphrev et al] (|2009l) used a parametric Bayesian 
analysis. The parametric methods used in the literature 
will generally underestimate the range of mass profiles 
consistent with the data. They also introduce systematic 
biases in the masses derived because of the assumptions on 

the profile shapes. 

W e use the sample of galaxies from IChurazov et al] 
(|2010l l and we discuss them and the properties of the hot 
X-ray gas they harbour briefly in Section [2] In Section [3] 
we describe the implementation of a new non-parametric 
Bayesian approach to obtain the total mass and circular ve- 
locity profiles from hydrostatic equilibrium. In Section [4] we 
show the results of tests of the new method and how to op- 
timise it. In Section [5] we show the total mass profiles and 
circular velocity curves we obtain by applying the method 
to the sample of galaxies, and the stellar and dark matter 
contributions we infer from published photometry and stel- 
lar population model mass-to-light ratios. In Section [6] we 
compare the individual circular velocity curves we obtain to 
previous determinations from X-rays and published dynam- 
ical models. We also examine the isothermality of the mass 
profiles of these galaxies, and look for any correlations that 
may exist between the stellar component, the total circular 
velocities and properties of the surrounding environment. 
We end with our conclusions in Section [7] 



2 THE SAMPLE 

We work on the sample of galaxies analysed in 
IChurazov et al] (|2010l ). Table [l] lists the galaxies in the sam- 
ple, the galaxy classification, the distances we assume, the X- 
ray luminosity and the velocity dispersion of the surround- 
ing environment. NGC 1399 is located at the centre of the 
Fornax cluster, NGC 1407 and NGC 5846 at the centre of 
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groups, and NGC 4472, NGC 4486 and NGC 4649 at the 
centre of sub-clumps in the irregular Virgo cluster. NGC 
4486 is at the ce n tre of the most massive sub-component. 
iDrinkwater et al.l (|2001 ) determined the velocity dispersion 
of the environment surrounding NGC 1399 from 92 galaxies 
within a radius of about 2.7° (930 Mpc), iTrentham et al.l 
ll2006l) from 35 galaxies w ithin about 900 kpc of NGC 1407 
and iMahdavi et al.1 (2005) from 87 galaxies within a radius 
of 1.8° (777 kpc) of NGC 5846. iGavazzi et al.1 l|l999h calcu- 
lated the velocity dispersion in sub-clumps containing NGC 
4472, NGC 4486 and NGC 4649 using 166, 62 and 58 galax- 
ies respectively. 

The sample galaxies harbour significant amounts of hot 
gas that except in the case of NGC 4486, appears relatively 
undisturbed, suggesting hydrostatic equilibrium. 



2.1 Density and temperature profiles of the hot 
gas 

The X-ray spectrum of the hot gas surrounding X-ray 
bright ellipticals consists of many emission lines and con- 
tinuous emission primarily from the mechanism of thermal 
bremsstrahlung. To obtain information on the temperature 
and density profiles of the gas, the spectra are fitted with 
models that make assumptions about the absorption along 
the line-of-sight and the metal abundance in the gas. 

We use the deprojected temperatu re and density pro- 
files obtained bv lChurazov et al.l (|2010l ) from Chandra and 
XMM-Newton observations, shown in Figure[T] The compos- 
ite profiles benefit in the central region from the high spa- 
tial resolution of Chandra and in the outer parts from the 
large field-of-view and collecting power of XMM-Newton. 
The XMM-Newton profiles are not more spatially extended 
than the Chandra profiles as one would normally expect, 
because several pointings were used with the Chandra in- 
strument. 

Churazov et al. first spherically deprojected the X-ray 
spectra us ing a non-parametric l east-squares algorithm de- 
scribed in IChurazov et al.1 (J2003). This finds the set of 3-D 
spectra in spherical shells that project to best-fit the ob- 
served spectra, assuming a power-law decline in the emis- 
sion outside the maximum observed radius. The deprojected 
spectra were then fit with the single-temperature APEC 
code in XSPEC, which resulted in a determination of the 
3-D temperature and density for each shell. The abundance 
of heavy elements in the model was fixed at 0.5 solar in all 
shells. This was assumed because for cool gas with approxi- 
mately solar abundance of heavy elements, the contributions 
of continuum and lines are difficult to disentangle unambigu- 
ously, resulting in an anti-correlation between the emission 
measure and abundance. 

Figure [l] shows generally a good agreement between the 
Chandra and XMM-Newton profiles in the region of over- 
lap. NGC 1399, NGC 1407, NGC 4649 and NGC 5846 have 
deprojected temperature profiles that start around 1.0 keV 
in the central region, dip slightly towards 0.5 keV and then 
slowly increase outwards to about 1.5 keV. The temperature 
profile of NGC 4472 increases steadily from about 0.7 keV to 
1.4 keV while the temperature profile of NGC 4486 increases 
much more steeply from about 1 keV to 3 keV with a small 
dip in between. The density profiles are much smoother than 
the temperature profiles. They appear to be very similar to 



each other with a linear decrease in a log-log scale and are 
therefore close to power laws. 



3 NON-PARAMETRIC RECONSTRUCTION 
OF THE MASS DISTRIBUTION 

If the gas is in hydrostatic equilibrium, we can obtain the 
mass distribution from the deprojected gas temperature and 
density profiles. Here we describe a new non-parametric 
Bayesian approach to obtain the most probable mass profiles 
within some confidence range. 



3.1 Hydrostatic equilibrium 

Hydrostatic equilibrium relates the 3-D temperature and 
density profiles of the gas to the 3-D mass distribution of the 
galaxy. If the gas is relatively undisturbed, we can assume it 
is in hydrostatic equilibrium. Balancing the pressure of the 
gas with the gravitational forces acting on it in a spherical 
system gives: 



d$(r) 1 dP(r) 



(It 



p(r) dr 



(1) 



where p is the gas mass density, P the gas pressure and r 
is the 3-D radius from the centre of the galaxy. We assume 
that the gas is ideal and therefore P = nksT and p = pm p n, 
where n is the particle number density of the gas, and p = 
0.61 is the average gas particle mass in terms of the proton 
mass, m p . This value of p corresponds to a helium number 
density of 7.92 x 10 -2 and 0.5 solar abundance of heavier 
elements. Now we have related the temperature and density 
of the gas to the gravitational potential it resides in. We can 
express Equation fT} in terms of the circular velocity curve, 
a distant-independent measure of the mass, V^ = rd^/dr. 
The circular velocity at some radius in a galaxy is the orbital 
velocity of a star at that radius on a circular orbit in the 
same gravitational field. Equation |T} becomes: 



y2 k b T rflnP 

prrip dlnr 



(2) 



It can be seen that to obtain the circular velocity curve, 
the temperature and gradient of the logarithmic pressure 
is required. This is more convenient than the usual form 
of this equation, where both the temperature and density 
derivatives are required. 



3.2 Bayesian approach 

We would like to find the most probable logarithmic pressure 
gradients (and then circular velocity curve from hydrostatic 
equilibrium) within some confidence range given the depro- 
jecte d temperature and density p rofiles measured from the 
data. iMerritt fc Tremblavl (1 19941 ') showed that it is best to 
treat these problems non-parametrically because paramet- 
ric methods are susceptible to systematic biases and under- 
estimate the confidence ranges. The procedure we develop 
can eventually incorporate the spectral fitting and depro- 
jection directly, but in a first step we concentrate here on 
generalising parametric methods on deprojected profiles to 
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Figure 1. Properties of the hot gas surrounding the sample of six X-ray b right elliptical galaxies: Deprojected temperature (upper 
plot) and density profiles (lower plot) were obtained bv lChurazov et al,l |2010) from the analysis of Chandra observations (red dots) and 
XMM-Newton observations (blue dots). The temperature is shown on a linear scale. The density and radii are shown on logarithmic 
scales. 



non-parametric ones. We adopt a Bayesian approach rather 
than the more conventional method of \ 2 minimisation for 
two reasons: Firstly, prior information such as the intrinsic 
smoothness of the profiles is easily incorporated. Secondly, 
the probabilistic nature of Bayesian methods means confi- 
dence ranges are more easily extracted. 

Let us describe the galaxy by some model M. Bayes' 
theorem tells us the probability of M given the deprojected 



temperature and density profiles X from the X-ray observa- 
tions: 



p(M | X) 



p(X 1 M)p(M) 



(3) 



Taking the logarithm of both sides gives: 



ln(p(M | X)) = \n(p(X | M)) + ln(p(M)) - ln(p(X)) (4) 



X-ray bright elliptical galaxies 5 



p(M | X) is the posterior probability of the model given 
the deprojected temperature and density profiles. p(X | M) 
is the likelihood probability C(M) of the deprojected tem- 
perature and density profiles given the model. p(M) is the 
prior probability of the model, indicating what we thought 
about the probability of that model before we knew anything 
about the deprojected temperature and density profiles. Fi- 
nally p(X) is the probability of those particular deprojected 
density and temperature profiles, but here only acts as a 
normalising factor because the profiles are fixed by the X- 
ray observations. To obtain the most probable model profile, 
we need to maximise the posterior probability p(M | X). In 
order to find all model profiles within some confidence range 
we have to find the shape of the posterior probability distri- 
bution. 

We define our model M = M(T, Ipg, P, V c ,m), where T 
is the temperature, Ipg is the logarithmic pressure gradient, 
P is the pressure, V c is the circular velocity and m is the en- 
closed mass. The pressure P is obtained by integrating the 
logarithmic pressure gradient Ipg. The circular velocity V c is 
obtained by applying hydrostatic equilibrium, and the mass 
m follows from V c 2 = Gm/r. The model T and P are com- 
pared with the deprojected Tdcpro and Pdepro obtained from 
the X-ray observations and deprojection. The model func- 
tions are defined on a grid of n mo d logarithmically-spaced 
radii r, finer than the grid of ndepro deprojected radii rdepro- 
We interpolate the functions to the grid of deprojected radii 
for the calculation of the \ 2 function below. We shall use 
j to denote an element on the model grid and k to denote 
an element on the deprojected grid. For example, T[j] de- 
notes the model temperature at the radius r[j] on the model 
grid and T[k] denotes the model temperature at the radius 
»"depro[fc] on the deprojected grid. 

We assume that the likelihood function C(M) takes the 
form of a multi-variate Gaussian: 



C(M) « exp(- X 2 /2) 



(5) 



X is the goodness-of-fit of the model temperature to the 
deprojected temperature and the model pressure to the de- 
projected pressure: 



model is very similar to the deprojected temperature and 
density at the deprojected radii but with varying values 
in between. Therefore applying smoothing constraints re- 
duces the model degeneracies by getting rid of unphysical 
solutions. We specifically penalise non-smooth solutions of 
the model temperature, logarithmic pressure gradient and 
resulting circular velocity curve because this is consistent 
within the hydrostatic framework. We define the smoothing 
prior as follows: 



p(M) oc exp(-S) 



(7) 



S is a weighted sum of the mean-square second derivatives of 
the temperature, logarithmic pressure gradient and circular 
velocity curve: 



si= £ r^ u ' " "\v ji '"^ u " ) o) 



S = c(Si +S 2 + S 3 ) (8) 
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-21nV c [j]+lnV c [j- 


1] 
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(11) 



h is the logarithmic interval between the model radii 
in an equally-spaced logarithmic radial grid and c = 
rc,de P roA/(n m od — 2) determines the weight of the smooth- 
ing term relative to the \ 2 term. It is defined so that for 
fixed A, the smoothing term is independent of ndepro and 
n mo d. A is specified by the user and Section 14.21 describes 
how we choose its optimal value. As the second derivatives 
are difficult to estimate at the boundaries, we omit their 
contribution to the smoothing prior. 

Combining the likelihood function and the priors we get 
the following expression for the logarithm of the posterior 
probability, which we call F: 



E 



PdcproM ~ P[k] 

e P [k) 



E 



r d0 p ro [fc]-r[fc] 

e T [k] 



(6) 
where ep and er are 1-cr statistical errors on the deprojected 
pressure and temperature profiles respectively. This choice 
for the likelihood function assumes a Gaussian error distri- 
bution and can be generalised if needed. To enforce physi- 
cally acceptable solutions, we impose the following boundary 
conditions as priors on the model temperature and logarith- 
mic pressure gradient profiles: 



(i) T > 

(ii) V c £3t: 



Ipg^O 



We also need to define a smoothing prior because un- 
like in a parameterised model, the best-fit non-parametric 
model would go exactly through all the deprojected tem- 
perature and density points resulting in noisy profiles. Also, 
as the model grid is finer than the deprojected grid, with- 
out smoothing there would be many solutions where the 



F = \n(p(M | A-)) = -i- -S + N 



(12) 



TV is a normalisation constant, which we set to zero for con- 
venience. 



3.3 The distribution of posterior probabilities 

Now that we have defined our model and the posterior prob- 
ability of a particular model given the X-ray information, 
we need an initial model and a procedure to generate new 
models in order to find the distribution of posterior proba- 
bilities. Here we will discuss our initial model and the 'pro- 
posal function', a function used to generate a new model. 
We then describe an algorithm we implement, optimised to 
find the most probable model. This is similar to the method 



of lMagorrianl (|1999T ). who found the most probable 3-D dis- 
tribution of stars given a projected surface-brightness dis- 
tribution. We also describe a second algorithm based on the 
Metropolis-Hastings sampling scheme in a Markov Chain 
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Monte Carlo (MCMC) approach, optimised to extract con- 
fidence ranges associated with a model. We define an 'iter- 
ation' as a change proposed by the proposal function that 
lies within the boundary conditions. 



3.3.1 The initial model 

The initial model is given by 

Minit(Tinit,/pfl , init,.Pinit, V c ,init , m ini t ) , calculated on the 
model grid. Ti n i t is determined by fitting a highly smoothed 
spline to the deprojected temperature Tdcpro- Ipginit is set 
to the slope of a straight-line fit between the logarithm 
of the deprojected pressure, Pdepro and the logarithm of 
the deprojected radii, rdcpro- V^mit and m^it follow from 
applying hydrostatic equilibrium (Equation ((2|), and the 
posterior probability of this initial model -Fi n it can be 
calculated using Equation (|12l) . 



a(i, i + 1) = min(l, r) where r — exp 



AF 
(AF) 



(14) 



A random number z is generated between and 1. If 2 ^ 
<x{i, i + 1), the proposed change is accepted and the disper- 
sion of the proposal function is decreased, <jt[J] = <jt[J]/u 
and o-ipg[j] — <Ji pg [j]/u, where u — 1.01 to allow for fine 
changes. Otherwise the proposed change is rejected and 
ar\j] = ctt\J] x u and oi Pg [j] = <Ji pg [j] x u. This update 
to the proposal function is made after each change to the 
model so that in general we follow a path that maximises 
the posterior probability, but still allow for the possibility of 
jumping away from a local maximum. The value for (AF) 
is updated every ni og steps. 

This procedure is continued until the model no longer 
significantly changes. 



3.3.2 The proposal function 

The proposal function q(Mi+i \ Mi) is a probability distri- 
bution that generates a model Mj+i from model Mi. We 
only make changes to the model T and Ipg and then cal- 
culate the resulting model P from integrating Ipg and the 
model V c and m by applying hydrostatic equilibrium. The 
proposal function is defined as: 



q(M, 



>+l 



Mi 



Pb}G T [j}G lpg [j} 



(13) 



P[j] = l/ n mod, is the probability of picking point j on the 
model grid. Gt[J] — p(Ti+i[j] \ Ti[j]) is the probability of 
going from Ti\j] to T i+1 [j] and Gi pg [j] = p(lpg i+1 [j] | lpgi[j]) 
is the probability of going from lpgt[j] to lpgi+i\j], at the 
chosen grid point j. Gt[J] and Gi pg [j] are described by Gaus- 
sians centred on Ti[j] and lpgi[j] respectively and with dis- 
persions o-r[j] and cri pg \j]. ot is initially set to the root- 
mean squared deviation between Tinit and T dGpro , and ai pg 
is initially set to the root-mean squared deviation between 
Ipginit and a two-point estimate of the logarithmic pressure 
gradient. 



3.3.3 Obtaining the most probable solution (non-Markov 
mode) 

To obtain the most probable model we want to get to the 
maximum of the posterior probability distribution as quickly 
as possible. We first relax the initial model to mitigate any 
bias introduced in choosing it. In this phase the proposal 
function is used to make a change to the model at a random 
point on the model grid. If the boundary conditions are met, 
then the proposed change is accepted and the new posterior 
probability F and change in posterior probability AF are 
calculated. After n rc i ax iterations, the relaxation phase ends 
and the average change in the posterior probability (AF) is 
calculated. In the relaxation phase, it is possible for {AF) 
to be either negative or positive as we accept all changes 
subject to the boundary conditions. 

In the improvement phase, a change to the model is 
proposed via the proposal function and if the boundary con- 
ditions are met then an acceptance probability a is specified 
for going from model Mi to model Mi+i: 



3.3.4 Obtaining the confidence range (Markov mode) 

To obtain the confidence range we need to probe the shape 
of the posterior probability distribution. The most efficient 
way to do this is using the Metropolis-Hastings algorithm, 
which generates models with a probability equal to the pos- 
terior probability. In our implementation of this algorithm 
we have three phases. The first is a relaxation phase of the 
initial model as in Section [3.3.31 with n rc i ax iterations. The 
second is a tuning of the proposal function, where the pro- 
posal function is updated as in the second phase in Section 
13.3.31 This phase ends after nt une steps, defined such that 
the acceptance rate in the third phase is 23%. This is recom- 
mended as the most efficient rate for probing the shape of the 
posterior probability distributions in high-dimension al im- 
plcmentations of the Metropolis-Hastings algorithm (jLiddlel 
12009(1 . 

In the third phase we want to map out the shape of the 
posterior probability distribution around the maximum. The 
proposed changes must therefore satisfy 'detailed balance': 



p(Mi\X)p(M^ 



Mi 



■ p(Mi+i \X)p(Mi\M i+1 ) (15) 



p(M i+ i I Mi) = q(M i+ i I Mi)a(i,i + 1) andp(M* j M i+1 ) = 
q(Mi I Mi+i) a(i + 1, i), where the proposal function q and 
acceptance probability a were defined in Equations (|13jl 
and (|14p respectively. This is ensured by no longer changing 
the proposal function and replacing the definition of r from 
Equation ([T4")l with: 



r = exp AF 



(16) 



r is called the Metropolis ratio. This will produce a sequence 
of models called a Markov chain. This is continued until the 
Markov chain converges to the posterior probability distri- 
bution. The initial ~ 25-50% of the accepted changes is 
called the 'burn-in phase' and is discarded because they do 
not reflect the posterior probability distribution. The con- 
vergence of the Markov chain can be checked by looking at 
the distribution of model values at each grid point and ensur- 
ing that subsequent iterations produce insignificant changes. 
Marginalisation with a Markov chain is trivial because 
the density of any grid point in the Markov chain will be 
proportional to the posterior distribution marginalised over 
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Figure 2. Application of the non-parametric analysis to a test model: The black lines show the (i) temperature, (ii) density and (iii) 
circular velocity curves for our test model. The filled black circles show the pseudo temperature and density profiles. The remaining lines 
show recovered profiles from pseudo deprojected temperature and density profiles generated from the test model, assuming A = 0.001 
(red, dotted), A = 0.5 (green, dashed) and A = 100 (blue, dash-dotted). 



the other variables. The shape of the marginal posterior den- 
sity at each grid point will tell us the median and confidence 
ranges, and the expectation value is simply a mean of all the 
values generated in the Markov chain. 



4 TESTS 

We need to test the ability of the two methods to recover 
the circular velocity curve of a galaxy and to choose values 
for the parameters n re i ax , niog, r&burnim ^iter and A that we 
have introduced in Section [3.21 13.3.31 and [373741 

For the relaxati on phase we set n r e i ax = 3(2n mo( j), sim- 
ilar to that used by iMagorrianl (J1999). In the non-Markov 
mo de we set n i og = 8(2n mo d), again similar to that used 
by IMagorrianl (1999). In the Markov mode we set Hburnin 
to 25% of the total number of accepted changes to ensure 
we are only sampling the posterior probability distribution. 
We choose the total number of accepted changes ni tcr as 
the number of iterations after which the model no longer 
appears to change significantly. 

To test the two modes of operation and calibrate the 
number of accepted changes mtm and the smoothing param- 
eter A, we define a model of a typical X-ray bright elliptical 
galaxy, from which we draw pseudo 3-D temperature and 
density profiles. 



4.1 The test model and pseudo temperature and 
density profiles 

The test model is defined by a 3-D temperature pro- 
file, 3-D density profile and a circular velocity curve 
[Ttesti "test, Vc.tost] typical for a X-ray bright elliptical galaxy. 
These profiles are shown by the black lines in Figure [21 For 
the circular velocity curve we use the dynamical model of 
NGC 5846 bv lKronawitter et all l|200Ch . We assume the tem- 
perature Ttcst is constant at 0.65 keV and then calculate 
the corresponding Ipgtest and density (n tcs t) profiles. n? eat is 
projected and pure Poissonian statistics typical for XMM- 
Newton is added . The projected densit y profiles are then 
deprojected as in lChurazov et al.l |2008|) and a deprojected 
density profile is obtained. We added Gaussian random de- 
viates to the temperature assuming a constant error of 0.03 
keV to complete the set of pseudo deprojected profiles shown 
with filled black circles in Figure [2] 




0.001 0.010 0.100 1.000 10.000 100.000 
A 

Figure 3. Residuals between the circular velocity curve of the 
test model and the circular velocity curves recovered by the non- 
parametric analysis for different values of the smoothing param- 
eter, A. 



4.2 Choice of optimal smoothing parameter 

We use our method in the non-Markov mode to obtain the 
best-fit circular velocity curve for different values of A from 
the pseudo deprojected temperature and density profiles. 
We find that niter = 10 6 iterations is sufficient to reach a 
best-fit model. 

A plot of the root-mean squared deviation R between 
the recovered circular velocity curve of the test model and 
the input circular velocity curve, against the smoothing pa- 
rameter A, is shown in Figure [3] where: 



R 



E(Vc,t 



Vc. 



^raod 



1/2 



(17) 



Figure [3] shows that R has a minimum at around X opt = 0.5. 
A lower smoothing leads to unphysical fluctuations in the 
recovered model and a higher smoothing results in a recov- 
ered model that does not fit the temperature and density 
profiles of the test model well. 

Figure [2] shows the resulting best-fit temperature, den- 
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Figure 4. The circular velocity values V c generated at each grid point j from our non-parametric analysis of the pseudo deprojected 
temperature and density profiles: The filled black circles show the binned number density (proportional to the marginal posterior density) 
of V c values at each grid point j for (a) niter = 10 4 and (b) ni tcr = 2 X 10 6 and the green lines show the best-fit Gaussians. 
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Figure 5. Application of the non-parametric analysis to the pseudo deprojected temperature and density profiles: Going counter-clockwise 
from the top left are the (i) logarithmic pressure, (ii) logarithmic pressure gradient, (iii) temperature, (iv) mass and (v) circular velocity 
profiles. The left three profiles are on a logarithmic radial scale, while the two profiles on the right are on a linear radial scale. The red 
lines correspond to the input test model, the green lines correspond to the recovered expected profiles and the blue lines correspond 
to the recovered median profiles. The grey region shows all generated points, the pink region shows the 95% confidence range and the 
orange region shows the 68% confidence range. 
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Galaxy 


Xt 


XlnP 


M 25 (IO^Mq) 


V c ,25 (km/s) 


c 


C<io 


C>io 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC 1399 


2.7 


10.1 


13.2 ±1.1 


475.0 ± 20.9 


1.15 ±0.03 


1.16 ±0.06 


1.60 ±0.06 


NGC 1407 


1.4 


3.8 


21.6 ±7.0 


608.9 ±97.9 


1.26 ±0.07 


0.88 ±0.05 


1.26 ±0.09 


NGC 4472 


2.7 


1.9 


14.1 ±1.0 


491.8 ±17.0 


1.19 ±0.02 


1.02 ±0.02 


1.59 ±0.01 


NGC 4486 


16.2 


28.9 


15.7 ±0.6 


520.3 ± 9.3 


1.36 ±0.01 


1.60 ±0.04 


1.37 ±0.01 


NGC 4649 


3.0 


9.4 


16.4 ±4.3 


530.4 ±69.7 


1.13 ±0.01 


1.08 ±0.01 


1.38 ±0.01 


NGC 5846 


2.3 


3.7 


12.5 ±1.9 


463.2 ± 35.0 


1.17 ±0.05 


0.88 ±0.05 


2.16 ±0.02 



Table 2. Derived properties of the sample of six X-ray bright elliptical galaxies: (1) Galaxy name, (2) \ 2 per data point between 
deprojected temperature and model temperature, (3) x 2 P er data point between deprojected logarithmic pressure and model logarithmic 
pressure, (4) enclosed mass at 25 kpc and associated 95% confidence range, (5) circular velocity at 25 kpc and associated 95% confidence 
range, (6) power-law index of mass profiles on a log-log scale and associated 1-cr error, (7) power-law index of mass profiles on a log-log 
scale fitting only until 10 kpc and associated 1-cr error and (8) power-law index of mass profiles on a log-log scale fitting only above 10 
kpc and associated 1-cr error. 



sity and circular velocity curves for the minimum (0.001), 
maximum (100) and optimal (0.5) values for A. 

4.3 Confidence range on recovered circular 
velocity curve 

To obtain the range of circular velocity curves within some 
confidence level, we have to ensure that our method is cor- 
rectly probing the shape of the posterior probability. Assum- 
ing the optimal smoothing parameter found above, we now 
use the Markov mode on the pseudo profiles. n tun c = 2000 
results in an acceptance rate of about 23% in the subse- 
quent proposed changes. We use a higher ruter = 2 x 10 6 , 
above which the profiles change insignificantly and choose 
"■bumin = 5 x 10 5 , above which the expectation values calcu- 
lated change insignificantly. In Figure|4]we have plotted how 
the circular velocity points at each radius are distributed in 
the Markov chain (i.e. after the burn-in phase), to illustrate 
the difference between using a ni tor that is too low and one 
that is about right. In the left plot where ni tor = 10 4 and 
"-bumin = 5000, the distribution of circular velocity points 
at each radius (i.e. the marginal posterior probability of the 
circular velocity curve at each grid point) have not settled to 
a smooth distribution, and therefore one would suspect that 
the chain is not yet properly sampling the posterior prob- 
ability distribution. The plot on the right however shows 
that running the Markov chain for longer results in smooth 
marginal posterior probability distributions that are well de- 
scribed by Gaussians (illustrated by the green lines). 

From Figure [4jb) the median value for the circular ve- 
locity at each radius is the circular velocity for which the 
marginal posterior probability is highest. The expectation 
value is the mean of all the circular velocity values chosen at 
that radius and the associated c% confidence range is given 
by the extrema of c% of the circular velocity values with the 
highest marginal posterior probabilities. Figure [5] shows the 
results of calculating these for all the model profiles. One 
can see that the median and expected profiles lie on top of 
each other, showing that the marginal posterior probability 
distributions at each radius are symmetrical. It can also be 
seen that the test model profiles lie within or on the bound- 
ary of the 68% confidence range. We repeat this for different 
sets of pseudo deprojected density and temperature profiles 
and we are convinced that the procedure is able to recover 
the circular velocity curve without a systematic bias. 



5 TOTAL, STELLAR AND DARK MATTER 
MASS PROFILES 

Now that we have tested the method, we can apply it to the 
sample of six X-ray bright elliptical galaxies to derive total 
mass and circular velocity profiles (Figure [f5] and Table [5}. 
We then estimate the stellar contribution to the total mass 
profiles using optical photometric data and stellar popula- 
tion models. Assuming negligible gas mass, we subtract the 
stellar mass from the total mass profiles to infer the dark 
matter mass fractions. 



5.1 Total mass profiles and circular velocity 
curves 

The fits to the deprojected logarithmic pressure and temper- 
ature profiles and the derived logarithmic pressure gradient, 
circular velocity and mass profiles are shown in Figure [6] for 
the sample of six galaxies. Columns (2)-(7) in Table [2] list 
the x 2 P er data point of the expected models, the enclosed 
mass and associated 95% confidence range at 25 kpc and 
the circular velocity and associated 95% confidence range at 
25 kpc. We choose 25 kpc because the deprojected profiles 
for all the galaxies extend at least as far as this radius. We 
quote the x 2 P er data point rather than the x 2 per degrees 
of freedom because the latter is unknown. The degrees of 
freedom is equal to the number of constraints subtracted by 
the number of parameters. In our case the number of con- 
straints is the number of data points plus the number of 
constraints introduced by the smoothing, which is difficult 
to define. 

Generally, the fits in the mean are very good. The x 2 
per data point between the deprojected temperature and 
model temperature ranges from 1.4 to 16.2. The x 2 P er 
data point between the deprojected logarithmic pressure 
and model logarithmic pressure ranges from 3.7 and 28.9. 
The x 2 values are not comparable to the usual reduced x 2 
which is ~ 1 at the 1-cr level, because the reduced \ is 
normalised by the number of degrees of freedom. Our x 2 
values may also be higher than expected because we have 
only considered statistical errors on the deprojected temper- 
ature and pressure profiles. There is also correlated/anti- 
correlated scatter in the deprojected temperature profiles 
betwee n adjacent points as a result of the deprojection pro- 
cedure (|Churazov et al.l l2008l. AGN activity driving shock 
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Figure 6. Application of the non-parametric analysis to the sample of six X-ray bright elliptical galaxies: Going from top down for each 
galaxy are shown the (i) logarithmic pressure, (ii) logarithmic pressure gradient, (iii) temperature, (iv) circular velocity and (v) mass 
profiles. The top three profiles are on a logarithmic radial scale, while the bottom two are on a linear radial scale. Black points on the 
logarithmic pressure and temperature plots show the deprojected profiles and associated statistical errors. The green and blue profiles 
(often indistinguishable) correspond to the expected and median profiles respectively. The grey region shows all generated points, the 
pink region shows the 95% confidence range and the orange region shows the 68% confidence range. The vertical black dotted lines on 
the circular velocity plots are at lR e , given in Table [3] 
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Galaxy 


Source 


Band 


Seeing (") 


R e (kpc) 


Source 


R2 (", kpc) 


M/L(M Q /L @ ) 


Source 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 1399 


soo 


B 


0.05 


4.17 


K10 


2.9, 0.28 


9.7 


KrOO 


NGC 1407 


Sp08 


B 


1 


3.32 


H09 


2.7, 0.38 


1.2 


Z07 


NGC 4472 


K09 


V 


0.05 


16.34 


H09 


3.5, 0.28 


6.0 


TOO 


NGC 4486 


K09 


V 


0.05 


31.44 


H09 


7.6, 0.59 


7.0 


G09 


NGC 4649 


K09 


V 


0.05 


10.52 


H09 


4.2, 0.34 


7.8 


TOO 


NGC 5846 


KrOO 


V 


0.2 


9.70 


KrOO 


3.3, 0.39 


7.2 


TOO 



Table 3. Photometric data for the sample of six X-ray bright elliptical galaxies: (1) Galaxy name, (2) source of surface-brightness 
profiles, (3) band of photometry, (4) seeing values , (5) effective radii, (6) source of effective radii, , (7) radius at which the optical 
surface-brightness profile has a gradient of 2, (8) stellar population mass-to-light ratios converted to the pho tometric band in c olumn 
(3) and distances fr om Table [T] and (9) source of stellar population mass-to-light ratios. References are SOO llSaglia et al.ll20001 ). Sp08 



Spolaor et al. 



Hopkins et al 



2008), 
2009) 



- , .. - - al. 21 

Z07 jZhang et alj|2007h . TOO JTrager et alj|2000h . G09 l|Gebhardt fc Thomadl2009h . 



K09 (Kormcndy et al. 



20091) . KrOO |K ronawi tter et a.l.ll2000f), K10 (J. Kormendy, private communication), H09 



waves into the ICM (e.g. iForman et al.ll2007f ) can also result 
in systematic undulations in the measured temperature and 
density profiles. The models do not fit the deviations caused 
by the correlated/anti-correlated errors and by systematic 
errors because we have calibrated the smoothing to a test 
model of an ideal massive elliptical galaxy. If the statistics 
of the data are low as in the case of NGC 1407, then the 
error bars cover the correlated/anti-correlated scatter and 
the % 2 values are low. The outstanding statistics (~ 0.5 Ms 
of Chandra observation) and very high surface brightness 
of NGC 4486 enables deprojected temperature and density 
profiles to be obtained in much narrower bins than for the 
other galaxies. This allows systematic deviations to be re- 
solved, but as the errors bars are small, they do not cover 
the correlated/anti-correlated scatter or the systematic de- 
viations. As a result the \ 2 values for NGC 4486 are on 
average about 5 times higher than for the other galaxies. 

The general shape of the derived circular velocity curves 
can be summarised by a steep rise in the centre (except in 
NGC 1407) followed by a slight dip at a few to 10 kpc and 
then a more gentle rise outwards, reaching circular velocities 
of 463-609 km/s at 25 kpc. The outward rise in the circular 
velocity curves is a combined result of temperature profiles 
that rise outwards and logarithmic pressure gradients that 
increase in magnitude outwards. The dips in the circular 
velocity curves are a result of dips in the logarithmic pressure 
gradient around the same radius. These dips may reflect true 
changes in the mass distribution or they may be a result of 
shocks in the gas propagating outwards. Shock waves can 
be produced by an unsteady outflow of rel ativistic plasma 
from a central black hole, as in NGC 4486 (jChurazov et al.l 
I2008T ). The shockfront is characterised by a sharp increase in 
the pressure inwards and then a decrease in the rarefaction 
region behind the shock. This manifests itself as a dip in the 
circular velocity curve at the position of the shockfront in 
NGC 4486. 

The enclosed mass profiles of all galaxies increase out- 
wards with an increasing gradient as expected from the 
increasing circular velocity curves, pointing towards non- 
isothermal mass profiles. The 95% confidence ranges in the 
circular velocities and enclosed masses at 25 kpc are reason- 
ably small and are smallest for the galaxies with the best 
statistics (NGC 4472 and NGC 4486) because the x 2 values 
are higher as discussed above. 



5.2 The stellar mass contribution 

To estimate the stellar contribution to the total circular ve- 
locity curve, we use surface-brightness profiles from various 
sources in the literature, given in column (2) in Table [3] and 
shown in Figure [7J a). We deproject the surface-brightness 
profiles assuming a spherical stellar distribution and a power 
law outside the range of the data. We then integrate the 
resulting 3-D luminosity density profile to obtain the lumi- 
nosity profile of the stars. Assuming constant stellar mass- 
to-light ratios with radius (columns (8) and (9) in Table [3} , 
we obtain stellar mass profiles and stellar circular velocity 
curves from the luminosity profiles. 

In order to find a scaling radius that best represents 
the luminosity of the stellar component, we examine two 
alternatives. The first is the mean effective radius, given 
in Table [3] for the sample of six galaxies. For NGC 1399 
we use the determination of J. Kor mendy (private com mu- 
nication) based on photometry in ICaon et al,l II 19941 ) and 
lLauer et al.l i|2007f ). We calculate the mean effective radii 
from the major-axis effective radii and ellipticities for NGC 
14 07, NGC 4472 NGC 4486 and NGC 4649 determined 
by iHopkins et al.l (120091) from ph o tome try in lLauer et al.l 
J2007IHNGC 13 9 9).lBender et al.l (| 19881 ) (NGC 1407) and 
IKormendv et all (120091') (NG C 4472, NGC 4486 and NGC 
4649). IHopkins et al.l ( 20091 ) determined double-Sersic fits 
to the profiles but calculated the effective radii for the total 
profiles. We do not use their determin ation for NGC 139 9 
because it is based on photometry from lLauer et al.l (|2007f) . 
which only probes the central region. For NGC 5846, we use 
the m ean de Vaucouleur's effective radius from lFaber et al.l 
(J1989T ) . The effective radii of massive elliptical galaxies are 
large and therefore require extended surface-brightness pro- 
files with excellent sky subtraction , which we believe is tru e 
in the case of the photometry of IKormendv et al.l (2009), 
but not in the older photometry where the sky subtraction 
may not be so accurate. Therefore as an alternative, we also 
calculate R2 (column (6) in Table [3} , the radius at which 
dfi/dlogR = 2, where fi is in magnitudes/arcsec 2 . Figure 
|7Ja) illustrates the similarity between the general shapes of 
the surface-brightness profiles when the radius is scaled by 
R2, and the surface brightness by its value at i?2- We ob- 
tain values of R2 ranging between 2.7-7.6", corresponding 
to 0.3-0.6 kpc, just outside typical seeing values of 1-2". 
Figure [7J a) shows that R2 is located around the radius that 
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Figure 7. Luminous properties of the six X-ray bright elliptical galaxies: (a) Surface-brightness profiles scaled by i?2 (the radius where 
dfi/dlogR = 2) and by the surface brightness fi(R2). (b) Kormendy relation between the B-band luminosity Lb(30R2) and R2 (red 
filled circles) and between the B-band luminosity Ls{Re) and R e (black filled circles). log(i?2) has been shifted by 1 to the right. The 
straight red and black lines show the respective best-fit straight lines. 



separates the core where the surface brightness decays gen- 
tly, from the region where it starts decaying more steeply, 
and therefore may be mostly reflecting the core properties. 
To determine which scaling radius should be used as 
representative of the galax y's luminosity, we plot an analogy 
to the Kormendy relation l|Kormendvlll977l ) in Figure [TJb). 
We plot the B-band luminosity of the sample galaxies at 
3O-R2 against R2 and the B-band luminosity at 1R S against 
R e and fit straight lines through both. The best-fit straight 
lines have slopes of 0.74 and 0.93 and Pearson product mo- 
ment correlation coefficient (PMCC) values of 0.70 and 0.94 
respectively, supporting correlations at about 88% and 99% 
levels of significance for a sample size of six. The correlation 
between the effective radius and the luminosity of the stellar 
component is much stronger, and therefore we will use this 
to represent the size of the stellar component. 



5.3 Scaled total circular velocity curves and dark 
matter mass fractions 

Figure HJa) shows the scaled total circular velocity profiles 
of the sample of X-ray bright galaxies from the X-rays along 
with the 95% confidence ranges with dashed lines. The radii 
are scaled by R e and the total circular velocities are scaled 
by their value at R e . In the central R e , all the profiles rise 
steeply then fall again before rising more gently from 0.5- 
lR e . The outer rise in the circular velocity curves noted in 
Section r5.ll is more easily compared between the galaxies in 
this scaled plot. The curve of NGC 5846 rises most steeply. 
The curves of NGC 1399, NGC 1407, NGC 4472, NGC 4486 
and NGC 4649 have a similar rise between l-2R e but the 
curves of NGC 1399 and NGC 1407 keep rising outwards, 
while the curve of NGC 4486 levels off outside 2R e . 

The mass in gas calculated from integrating the gas den- 
sity profiles is < 0.04% of the mass in stars in all galaxies 
except in NGC 4486 where it is ~ 0.4%. Therefore we treat 
the gas mass as negligible and assume that the baryonic 
mass component in the galaxy is solely due to stars. Sub- 



tracting the stellar mass determined in Section T5. 21 from the 
total mass obtained from the X-rays gives the dark matter 
mass, shown as a fraction of the total mass in Figure [8jb), 
outside 0.6-R e . Below this radius the dark matter fractions 
we obtain for some of the galaxies are below zero, which 
could be a result of either an overestimate in the stellar 
mass-to-ratios or an underestimate in the total circular ve- 
locity curves derived from the X-rays, the latter of which 
we discuss in more detail in Section 16.21 The dark matter 
fractions range between 50-80% at 2R e ~ 6.5-59 kpc in all 
galaxies except NGC 1407 where it is ~ 35%. Even further 
out, the profiles appear to converge to a value between 80- 
90%, pointing uniformly to a massive dark matter halo. The 
larger range in the dark matter mass fractions further in is 
indicative of dark matter mass profiles that have different 
shapes. For example the dark matter mass profile is steep- 
est for NGC 4486, implying a more compact dark matter 
component. 



6 DISCUSSION 

In this section we compare the circular velocity curves we 
obtain from our analysis to those obtained from previous 
X-ray determinations and published dynamical models. We 
then address the issue of the isothermality of the circular 
velocity curves and how this compares with the literature. 
Finally we look at how the total circular velocity curves may 
be correlated with properties of the stellar components and 
the velocity dispersions of the environments to which the 
galaxies belong. 



6.1 Comparison with previous X-ray 
determinations 

Figure [5] shows the total circular velocity curves for our 
sample of galaxies along with their 95% confidence ranges 
in black solid and dashed lines. For comparison we have 
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Figure 8. Total and dark matter mass properties of the six X-ray bright elliptical galaxies: (a) Total circular velocity profiles scaled 
by R e in radius and the total circular velocity V c (R e ). (b) Dark matter mass fractions scaled by R e in radius. The black dotted line 
highlights the dark matter mass fractions at 2R e . 



overplotted various X-ra y determinations in the literature. 
The green lines are from iFukazawa et al.l (|2006T ) , who used 
Chandra observati ons of NGC 1399 and NGC 4472. The 
blue lines are from iNagino fc Matsushital (120091 ) who com- 
bined Chandra and XMM-Newton observations of NGC 
1399, NGC 4472 an d NGC 5846, and the pink line is from 
I Zhang et all (|2007h who combined Chandr a and ROSAT 
data i n NGC 1407. IHurnphrev et al.l (|2006T l (solid orange 
lines) , Humph rey et al.l l|2008l ) (dashed orange line) and 
IHurnphrev et al.l (|2009) (dotted orange line) used Chandra 
observations of NGC 1407, NGC 4472 and NGC 464 9 and 
finally, the cyan line is from iMatsushita et al.l ((2002) , who 
used XMM-Newton observations of NGC 4486. Density and 
temperature profiles of the hot gas were obtained from the 
observations and parametrised. The total circular velocity 
curves were then obtained by differentiating the profiles and 
applying hydrostatic equlibrium. 

Discrepancies in the outer slopes are most likely a result 
of differences in the spatial extent of the data used, as the 
last few points anchor the outer slope. Smaller-scale differ- 
ences are most likely to arise from differences in the methods 
used to obtain the mass profiles. The methods employed in 
the literature are all parametric and therefore the shapes 
of the mass profile will have a systematic bias. In NGC 
1407 however, the discrep ancy between the o uter s lope of the 
circular velocity curve of IHurnphrev et al.l |2006) an d that 
of our determination and the determination of IZhang et al.l 
(2007) is more signficant than between the other profiles. 
The source of this is unclear. 



6.2 Comparison with dynamical models 

Figure [U] also shows dynamical models and their confidence 
ranges from the literature, in red solid and dashed lines. 
The vertical red dotted line demarcates the radial extent of 
the data used in the dynamical modelling and we plot the 
dynamical models to twice this radius, except in the case of 
NGC 1407 where we plot the models until the radius deemed 



as reliable by the author. For NGC 139 9 and NGC 5846 
we use the spherical dynamical models of lKronawitter et al.l 
(2000). They built models by superposing spherical distri- 
bution functions and adjusting the weights of their contribu- 
tion to the model by comparing the projection of the model 
to measured surface-brightness profiles and long-slit kine- 
matics extending to 97" and 99" respectively. Although this 
corresponds to only 5-10 kpc, the models provide a useful 
comparison in the central regions. 95% confidence ranges are 
provided on the mod els. For NGC 1407, we com pare with the 
dynamical model of lRomanowskv et al.l (|2009T ). who solved 
spherical Jeans equations for a varying anisotropy profile 
and fit to a surface-brightness profile and globular cluster 
kinematics extending out to ~ 590". They provide a 68% 
confidence range on their model. For NGC 4472, NGC 4486 
and NGC 4649 we compare with the models of Shen & Geb- 
hardt et al. (2010, in preparation), Gebhard t et al. (2010, 
in preparation) and lShen fc Gebhardt! (J2010T ). who used ax- 
isymmetric Schwarzschild models to fit surface-brightness 
profiles, integrated kinematics and globular cluster data ex- 
tending to 1200", 500" and 533" respectively. 68% confi- 
dence ranges are provided on the models. 

Comparing the X-ray circular velocity curves of this pa- 
per to those from dynamical models, we find that outside the 
central kpc in NGC 1399, the X-ray circular velocity curve 
is up to 16% lower than the dynamical modelling circular 
velocity curve. Outside a radius of about 7 kpc though, the 
X-ray determination lies within the 95% confidence range 
on the dynamical modelling determination. In the region of 
overlap between the X-ray and dynamical circular velocity 
curves of NGC 1407, the X-ray circular velocity curve rises 
from about 500 to 600 km/s, while the dynamical circular ve- 
locity curve rises from about 270 to 320 km/s, showing that 
the X-ray velocity curve is almost twice as high throughout 
and therefore rises twice as steeply. The X-ray circular veloc- 
ity curve of NGC 4472 is up to 27% lower than the dynamical 
modelling circular velocity curve within about 12 kpc, but 
then rises to up to 20% higher outside this radius. Within 5 
kpc, the X-ray circular velocity curve of NGC 4486 is up to 
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Figure 9. Total circular velocity curves determined from X-rays, compared with those from previous X-ray determinations and published 
dynamical models: The X-ray velocity curves and associated 95% c onfidence ra nges determined in this pape r are shown by the solid and 
dashe d black lines respectively. Green lines l|Fukazawa_et L _alJ 200a) , blue lines I INaeino fc Matsush ita 2009), the pink line \ Zhang et al.l 
2007), orange solid, dashed and dotted lines llHumphrev et al] [2006. 2008, 2009), and the cyan line (Matsushita et al. 2002J) show circular 
velocity curves determined from X-rays in the literature. The best-fit dynamical models along with their confidence range (for more 
details see text) are shown with solid and dashed red lines. The vertical red dotted line shows the radial extent of the data used in the 
dynamical modelling. 



45% lower than the dynamical circular velocity curve, and 
outside this region, the X-ray circular velocity curve is up 
to 21%. For NGC 4649, the confidence ranges overlap from 
7 kpc outwards and the circular velocity curves agree very 
well outside 10 kpc. Within this radius the X-ray circular 
velocity curve is at worst about 21% lower than that of the 
dynamical model. In NGC 5846, the X-ray circular veloc- 
ity curve is up to 18% lower in the central 12 kpc than the 
dynamical circular velocity curve but then rises much more 
steeply compared to the dynamical circular velocity curve. 



6.2.1 Uncertainties in the X-ray analysis 

Here we discuss uncertainties in the X-ray analysis, which 
could explain the discrepancies we observe between the dy- 
namical and X-ray circular ve locity curves. 

One subject discussed in IChurazov et al.l i|2010l ) is the 
contribution of non-thermal components in the gas to the to- 
tal pressure. The existence of these components would lead 
to a systematic bias in our determination of the mass pro- 
file derived from hydrostatic equili brium as we have only 
considered the thermal gas pressure. IChurazov et al.l ([2008) 
listed several non-thermal components that are consistent 
with AGN-heated cores of X-ray bright elliptical galaxies ex- 
hibiting properties of a cooling flow region: (i) Cosmic rays 
and magnetic fields uniformly mixed with the thermal gas, 
(ii) cosmic rays and magnetic fields forming bubbles that are 
free of thermal gas and (iii) microturbulence in the thermal 



gas. If the fraction of the total pressure that is due to thermal 
gas ft is constant with radius, then cf> t = ft<t>true, where 4> t is 
the potential calculated from the thermal pressure and 0true 
is the true potential. This translates to ft = K? t /V^. 2 truc . As 
^ ft ^ 1) non-thermal contributions to the pressure are a 
possible mechanism for lowering the circular velocity curve 
calculated from the X-rays and therefore could explain the 
discrepancies we see in NGC 1399, NGC 4472, NGC 4649 
and NGC 5846 in the central ~ 10 kpc and throughout NGC 
4486. Calculating ft (for r > 1 kpc) where the X-ray circular 
velocity curve is lowest compared to the dynamical circular 
velocity curve gives upper limits on the local contribution of 
non-thermal components to the total pressure of 30%, 45%, 
69%, 32% and 33% in NGC 1399, NGC 4472, NGC 4649 and 
NGC 584 6. These local, inner va lues are higher than those 
quoted in IChurazov et al.l (|201C|) , where the ft is a global 
estimate from fits to the total potential. 

Another possible source of discrepancy could be in the 
assumptions made about the abundance profile in the deter- 
mination of deproject ed temperature and density profiles. 
IChurazov et al.l (|2010h assumed a flat abundance profile for 
their analysis but calculated the average effect on the cir- 
cular velocity curves if this assumption was relaxed. They 
found that the circular velocity values generally increase for 
their sampl e by a value of ~ 2% . This small effect is in agree- 
ment with IWerner et al.l (|2009l ), who combined very high- 
resolution XMM Newton data with Chandra data in NGC 
4636 and concluded that allowing for a varying abundance 
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profile does not change the deriv ed deprojected t e mper ature 
and density profiles signficantly. IChurazov et al.l i|2010f) also 
considered the effect of the unresolved population of low- 
mass X-ray binaries and estimated that this generally lowers 
the circular velocity curve by ~ 2%, which to some extent 
counteracts the effect of allowing an abundance gradient. 

We also consider the effect of the extrapolation in the 
deprojection of the spectra on the derived circular veloc- 
ity profiles. Th e depr ojected density profiles presented in 
IChurazov et al.l (|2010i ) and in Figure [T] in this paper as- 
sume a power-law extrapolation based on a power-law fit to 
the data. We consider two alternative extrapolations of the 
Chandra density profile of NGC 4472, to investigate whether 
we can explain the discrepancy between the X-ray and dy- 
namical circular velocity curves outside ~ 11 kpc. If the 
index of the power law used in the extrapolation of the den- 
sity profiles becomes less negative by 0.1 (i.e. the density 
profile becomes more shallow), the outermost Chandra cir- 
cular velocity point decreases by about 10 km/s, and effects 
further in are not discernable. Therefore this is not a realis- 
tic explanation for the discrepancy we observe. If instead the 
power-law index becomes more negative by 0.1 (i.e. density 
profile drops off more steeply) , then the effect is to increase 
the circular velocity of the last two Chandra circular velocity 
points by < 10 km/s. The effects further in are negligible. 
This would not be sufficient to explain why the X-ray circu- 
lar velocity curve is so much higher in NGC 1407 compared 
to the dynamical circular velocity curve. 

Shocks in the gas result in a local violation of hydro- 
static equilibrium and in spatial correlations in the pres- 
sure and circula r veloc ity curves fSection l5.1|l . In NGC 4486, 
IChurazov et al.l i|2008i ) showed a simulated shock at 2.7 ar- 
cmin (~ 13 kpc), which could correspond to the acute bump 
seen in the pressure and temperature profiles around this ra- 
dius in Figure [5] They found that such a shock leads to a 
dip in the potential of NGC 4486. We observe a similar dip 
in the total circular velocity of NGC 4486 close to this ra- 
dius. In the other galaxies, shockfronts are not obvious in 
the deprojected profiles, possibly due to the lower statistics 
compared to NGC 4486. In general, small-scale deviations 
from hydrostatic equilibrium are simply not fit due to the 
smoothing prior in our method, while a larger-scale violation 
results in a systematic effect. 

IChurazov et al.l (|2010i ) also discussed the effects of as- 
suming spherical symmetry in their determination of the 
gravitational potential from X-rays and found that on aver- 
age the effect is unlikely to exceed 7%. As the effect is sys- 
tematic and will either increase/decrease the average slope 
of the potential, the resulting change in the circular velocity 
is of the order of 4%. 

In earlier work. lChurazov et al.l (|2008l ) estimated biases 
introduced when fitting multi-temperature plasma with a 
single-temperature model by comparing profiles obtained 
in different energy bands. In the case of NGC 4486, they 
compared temperature and density profiles calculated from 
broad-band spectra to that from hard-band spectra and 
found that differences are small outside 0.4 arcmin, which is 
a couple of kpc. Within this region however, AGN activity 
can result in hot and cool gas structures superposed along 
the line of sight. NGC 4486 is the most disturbed galaxy in 
our sample but this source of uncertainty could also explain 



some of the discrepancies we find in the central few kpc in 
other galaxies in the sample. 

6.2.2 Uncertainties in the dynamical modelling 

Here we examine uncertainties associated with the dynami- 
cal modelling methods, which could be the source for some 
of the discrepancies we observe between the dynamical and 

X-ray circular velocity curves. 

In the dynamical models of iKronawitter et al.l (2000), 
the data used extends only to 5-10 kpc. However in the 
X-ray bright elliptical galaxies in our sample, the density 
profiles of the stars are shallow in the outer parts, which im- 
plies that there is a lot of stellar mass at large radii where 
the orbital structure remains poorly constrained. If addi- 
tionally the orbital structure of these o uter stars were very 



radial, as in the numerical simulations of Abadi et al 



2006) 



2009) 



and the dynamical models of Ide Lorenzi et al.l (2008 
(/3 = 1 — af/Vr ~ 0.5), there would be many stars on highly 
elliptical orbits traversing large distances in the galaxy. In 
the central regions they turn around and contribute to the 
LOS velocity dispersion more than compared to if they were 
on a circular orbit there. Therefore a higher mass profile 
could be inferred in the central region, compared to what 
would be inferred if the outer kinematics are also available. 
To estimate how much extra mass may be inferred in the cen- 
tral region we look at the axisy mmetric Schwarzschil d mod- 
els of Coma elliptical galaxies in lThomas et al.l 1)20071 ). They 
used kinematic data extending to betwe en 13-30 kpc, fur- 
ther o ut than for the dynamical models in lKronawitter et al.l 
(2000), and provide 68% confidence ranges on their circular 
velocity profiles. Without the outer kinematics, they would 
have inferred a less massive dark matter halo, and therefore a 
higher central stellar mass-to-light ratio. Let us assume that 
this is approximately given by the upper boundary of the 
confidence range on the mass profiles in the central region. 
Therefore the inner mass profile would be overestimated by 
about 50 km/s without the outer kinematic data, which is 
about 10% of the circular velocities we obtain from the X- 
rays. Therefore it could partly but not wholly explain the 
discrepanicies in regions where the dynamical circular veloc- 
ity curves are higher than the X-ray circular velocity curves. 

A similar argument is relevant in cases where the dy- 
namical models do not have central kinematic constraints 
on the orbital structure, such as in NGC 1407, where only 
globular cluster kinematics were used from 15 kpc outwards. 
Central kinematic constraints are invaluable as they are usu- 
ally of much higher quality than kinematic data further out. 
In radially anisotropic systems in particular, they not only 
place stringent constraints on the central dynamical struc- 
ture, but also constrain the orbital structure at least out to 
a few times their maximum projected radius. 

An assumption on the symmetry of the stellar distribu- 
tion in the dynamical models can also lead to a systematic 
bias. Let us consider the case where the stellar distribution is 
assumed to be spherical but in reality is prolate axisymmet- 
ric along the line-of-sight (LOS). Let us position the LOS 
along the z-axis, which will then be aligned with the long 
axis of the prolate ellipsoid. Therefore the velocity dispersion 
along the LOS is higher than in the plane of the sky. As we 
are measuring only the LOS velocity dispersions, a spher- 
ical dynamical model would overestimate the average 3-D 



16 P. Das et al. 



velocity dispersion (0-30) and therefore also the mass and 
circular velocity. For an isotropic velocity dispersion tensor, 
cr3D,sp ~ \/3crzz- If the stellar mass distribution is in fact pro- 
late along the LOS, then <T3D, P r = (2o"^ + er 2 2 )i , where x 
and y are along the short axes and a X x = a vv because of ax- 
isymmetry. The Tensor Virial Theorem relates the intrinsic 
axial ratios of the stellar mass distribution to its rotation and 
the anisotropy of the velocity dispersion tensor. Assuming no 
rotation in thi s galaxy gives <T3D, P r ~ <r zz (l + 2<j°' 9 )2 (e.g. 
lGerhardlll994T ). where q is the intrinsic flattening. There- 
fore K 2 true/V c 2 » a| D , pr /a 3 2 D , sh « (1 + 2g a9 )/3. This es- 
timate does not take into account that dynamical mod- 
els created assuming different shapes for the galaxies may 
prefer different orbit distributions, which influence the \ 2 
minimisation usually employed. Nevertheless, it is a use- 
ful way for us to determine approximately how prolate a 
galaxy has to be along the LOS to explain the discrepan- 
cies we observe. Looking at the radii where the dynami- 
cal and X-ray circular velocity curves of NGC 1399, NGC 
4472, NGC 4486, NGC 4649 and NGC 5846 are most dis- 
crepant, we find that for r > 1 kpc, maximum flattenings 
of q = 0.52, 0.29, 0.14, 0.48. 0.47 r espectively would be re- 
quired. [TrembIay_^Merritt] (l996) found that the distribu- 
tion of intrinsic, short-to-long axis ratios peaks at around 
0.75 for bright elliptical galaxies. Therefore it is unlikely 
that prolate axisymmetry along the LOS is the cause for 
the observed discrepancies. 

Finally, the mass profiles explored in the dynamical 
models are paramerised so that that stellar mass contribu- 
tion follows the shape of the stellar luminosity profile (i.e. 
the mass-to- light ratio is constant with radius). The dark 
matter mass contribution is usually given by non-singular 
isothermal spheres or NFW profiles. This results in a bias in 
the shape of the mass profiles and circular velocity curves, 
especially noticeable in the central regions where the shape 
of the circular velocity profile changes most. 

Based on our discussion on the uncertainties inherent 
in the derivation of circular velocity curves from X-rays and 
from dynamical models, we conclude that the regions where 
the X-ray circular velocity curves are lower compared to 
dynamical circular velocity curves in NGC 1399, NGC 4472 
and NGC 5846 can probably be attributed to a combination 
of: i) A contribution of non-thermal components to the pres- 
sure in the X-rays, ii) multiple-temperature components in 
the gas, iii) an incomplete spatial converage in the kinematic 
data used in the dynamical models, and iv) insufficiently 
general mass profiles in the dynamical modelling. As the dy- 
namical models for NGC 4486 and NGC 4649 incorporate 
extended kinematic data, the probable explanation is: i) the 
existence of non-thermal pressure components, ii) multiple- 
temperature components, and iii) insufficiently general mass 
profiles in the dynamical models only. The discrepancy be- 
tween the X-ray and dynamical circular velocity curves in 
NGC 1407 remains unclear. 



6.3 How non-isothermal are the mass profiles? 

In the previous section we carried out an individual compari- 
son between the circular velocity curves we obtained to those 
obtained in the literature. Here we examine the general prop- 
erty of rising circular velocity curves noted in Section 15.11 
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Figure 10. The non-isothermality of the mass profiles of the 
sample of six X-ray bright elliptical galaxies: Solid lines show 
the mass distributions in a double-logarithmic scale. Dotted lines 
show straight lines fitted to the mass distributions in this scale. 



To quantify how non-isothernral the circular velocity curves 
are, we fit a p ower-law model M oc r^ to the mass pro- 
files a s done bv lFukazawa et al.l (J2006j) and Humphrey et al.l 
(2009) from their analysis of XMM-New ton and Chandra 
observations, and lChurazov et al] (J2010j) on the sample we 
work with. Our determinations of £ are listed in columns 
(6)-(8), calculated by fitting to the total mass profile, the 
mass profile within 10 kpc and the mass profile outside 10 
kpc. Figure [10] shows the mass profiles in a log-log scale with 
the best-fit models fit over the whole radial range. In this 
plot, the slope of the best-fit lines is equal to the power- 
law index £. For a galaxy with an isothermal mass profile, 
we would expect f = 1. Fitting to the the total mass pro- 
file we find 1.13 < C < 1-36 and a mean (£) = 1.21. As 
the 1-cr errors range between 0.01 and 0.07, the sample of 
these six elliptical galaxies appear to have signficantly non- 
isothermal mass profiles. Fitting the mass profiles within 10 
kpc, we find 0.88 < C < L60 and a mean {() = 1.10. For 
NGC 1399, C does not change, for NGC 4486, < becomes 
much larger, possibly due to the shockfront at this radius, 
but for the remaining galaxies £ decreases to values between 
0.88 and 1.08. Calculating £ from the mass profiles outside 
10 kpc, we find higher values ranging from 1.26 < C < 2.16 
and a mean {£) = 1.56. This implies that the central 10 kpc 
of these galaxies may be isothermal, but outside this radius 
the circular velocity curves are rising. The rising circular ve- 
locity curves could be showing that we have probed the mass 
distribution sufficiently far out to observe the effects of the 
massive group/cluster-sized haloes in which these galaxies 
reside. 

iFukazawa et al] (J2006J) obtained mass profiles by apply- 
ing hydrostatic equilibrium to parametrised Chandra and 
XMM-Newton temperature and density profiles for a sam- 
ple of 53 elliptical galaxies. They fit power-law mass models 
outside 10 kpc, finding on average £ = 1.33 ± 0.33, which 
over laps with our de t ermin ation of 1.56. 

IChurazov et al] i|2010r ) found a mean (£) = 1.11 for 
NGC 1399, NGC 4472, NGC 4486, NGC 4649 and NGC 
5846, fitting their potentials between 0.1' and 5'. Over the 
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same radial range for the same galaxies, we obtain (£) = 
1.17, which is only ~ 5% higher. 

iHumphrev fc Buotd (|2010r ) derived mass profiles from 
a parametric Bayesian analysis of Chandra observations for 
a sample of 10 galaxies, groups and clusters. They obtained 
0.95 < C, < 1.8 fitting to the mass profiles between 0.2- 
10-Re but for the four galaxies in their sample, they obtained 
0.95 < C < 1-09, i.e., close to isothermal. We believe that the 
differences between their results and ours is due to: i) Calcu- 
lation of £ over different radial ranges from ours and ii) envi- 
ronmental effects, i.e. our galaxies are located specifically at 
the centre of sub-groups, groups or clusters while their galax- 
ies, though belonging to similar environments, are generally 
not at the centre. The radial range is important in the case of 
NGC 4649, the only galaxy common to both samples. They 
used data until about ~ 17 kpc, converted to our assumed 
distance. If we look at Figure[T]we can see that there is an ex- 
tra deprojected temperature and density point that we have 
from the XMM-Newton analysis, which increases the gradi- 
ent of the mass profile. There is also reason to believe this 
increase because of h ow similar our derived cir cular velocity 
curve is to that from lShen fc Gcbhard tl 1120101) discus sed in 
Section H] For NGC 1332 IHumphrev fc Buotd IJ2010h used 
data only until ~ 20 kpc, which also may not be sufficiently 
extended (our profiles all extend to at least ~ 25 kpc and 
in the case of NGC 4486, until almost 110 kpc.) NGC 720, 
NGC 4261 (for which they have data until ~ 31 and ~ 34 
kpc) and NGC 1332 are not located at the centre of their re- 
spective environments. Therefore even if the profiles were as 
radially extended as ours, we would not expect the circular 
velocity curves to rise as much as for our sample galaxies, 
which are all located at the centre of th eir environments. For 
the gr oups and clusters in the sample of lHumphrev fc Buotd 
(J2010J), they found mass profiles that increase more steeply 
than isothermal mass profiles, therefore supporting the hy- 
pothesis that we are partially probing the region where the 
group/cluster potentials start dominating. The comparison 
between calculating £ within 10 kpc with calculating it over 
the whole radial range of the deprojected temperature and 
density profiles also supports this hypothesis. 

Gerhard et alj u200lj) analysed the dynamical models of 

iKronawitter et al.l (|2000r t and found that their circular ve- 
locity curves are flat to within 10% for R > 0.2R e . These 
models use photometric and kinematic data extending to 
less than lR e for 9/21 of the galaxies, between l-2_R e for 
6/21 of the galaxies a nd 2-3R e for 6/21 of the galaxies. 
iKoopmans et al.l 1120091 ) solved constant-anisotropy spheri- 
cal Jeans equations for 58 strong-lens early type galaxies us- 
ing stellar velocity dispersion constraints and lensing-based 
total masses within the central effective radius. They found 
an average logarithmic mass-density slope of -2.085 (for an 
isothermal mass profile, one would expect -2.0) and that the 
dependence of this result on the value of the anisotropy pa- 
rameter /3 is small. Both studies used elliptical galaxies that 
cover a range of luminosities and are not constrained to re- 
side at the centre of sub-group/group/cluster environments. 
This means that for many of their galaxies, we have no rea- 
son to expect a rising circular velocity curve and for those 
that we do, the mass profiles may not be probing far out 
enough and/or the lack of extended photometric and kine- 
matic constraints could be biassing the inner mass profile as 
discussed in Section ^, 2. 21 
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Figure 11. Correlations between the stellar component, the total 
circular velocity, and the environment: (a) Tully-Fisher relation 
between the luminosity of the stellar component at lR e and the 
total circular velocity from the X-rays at lR e (black filled cir- 
cles) and associated 95% confidence range (horizontal bars). Also 
shown are the results of considering the average effect of the un- 
resolved population of low-mass X-ray binaries and an abundance 
gradient (black open circles), (b) Total circular velocity from X- 
rays at lR e and the velocity dispersion of the environment (black 
filled circles), (c) Luminosity of the stellar component at lR e 
against the velocity dispersion of the environment (black filled 
circles). Solid black and dashed black lines show straight-line fits 
to the black filled circles and open circles respectively. 
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6.4 Correlations between the total circular 

velocity curves, the stellar component and the 
environment 



The Tully-Fisher relation links the luminous component to 
the total mass of the galaxies and for spiral galaxies is usu- 
ally shown between the total luminosity and the flat part 
of the circular velocity curve. We plot an analogous relation 
for our sample between the total B-band luminosity at R e 
against the total circular velocity from the X-rays at R e , 
shown with black filled circles in Figure 1111 We have also 
plotted with black open circles, the approximate total cir- 
cular velocities at these radii if we take into account the 
average effects of the unresolved population of low-mass X- 
ra y binaries and varying abundance profiles, as estimated 
bv lChurazov et al.l (J2010I ). The solid black and dashed black 
lines show straight lines fit to the black filled circles and 
black open circles respectively, considering the errors in the 
circular velocity. The lines have slopes of 3.77 ± 0.03 and 
4.01 ± 0.03, very similar to the value of 4 in the Tully-Fisher 
relation for spiral galaxies and to the value of ~ 4 (private 
communication with J. Thomas) obtained fo r the sample of 
Coma cluster elliptical galaxies analysed in iThomas et al.l 
(2009). PMCC values for the fits are 0.86 and 0.85, there- 
fore supporting a linear correlation at a ~ 95% level of sig- 
nificance for a sample size of six. This implies that systems 
with more luminous stellar components tend to have higher 
circular velocities. 



Figure lllf b) shows the relation between the circular 
velocity at R e and the velocity dispersion of the surround- 
ing environment from column (6) in Table [T] The solid line 
shows the best-fit straight line with a slope of 0.53±0.02 and 
a PMCC value of 0.92, supporting a linear correlation at a 
98% level of significance. This strongly implies that central 
galaxies residing in hotter environments (as measured by the 
velocity dispersions) have higher circular velocities. Figure 
lllf c) shows the logarithm of the luminosity at R e against 
the logarithm of the velocity dispersion of the surrounding 
environment. The best-fit straight line has a PMCC value 
of 0.83, supporting a linear correlation at a level of ~ 95% 
signficance. This implies that the luminosity of the galaxy 
residing at the centre of a sub-clump, group or cluster en- 
vironment, is related to the velocity dispersion of the envi- 
ronment. 



In environments where the local velocity dispersions are 
higher, the densities (~ a 2 /r 2 ) are also greater, as the radii 
(0.7-1 Mpc) over which these velocity dispersions were cal- 
culated are comparable. Therefore one would expect that a 
larger number of systems would fall onto the central galaxy, 
resulting in a more lumino us stellar component and a more 
massive dark matter halo. iBrough et al.l (|2007r i found that 
more X-ray luminous systems are intrinsically brighter in 
the K-b&nd compared to less luminous X-ray luminous sys- 
tems for their sample of three brightest group galaxies and 
three brightest cluster galaxies. In denser environments, one 
would expect the gas to be hotter and denser, and therefore 
more X-ray luminou s, thereby corroborating the results of 
iBroughet ai] | |2007l ). 



7 CONCLUSIONS 

In this paper, we describe a new non-parametric Bayesian 
approach to obtain mass distributions and associated con- 
fidence ranges from temperature and density profiles of hot 
gas in hydrostatic equilibrium. The method is able to suc- 
cessfully reconstruct the mass distribution of a test galaxy 
just within a 68% confidence range of the recovered models. 
We assume a smoothing prior to ensure unique and physical 
mass distributions and calibrate this on the test galaxy. 

We apply the procedure to the sample o f six X-ray 
bright elliptical galaxies from IChurazov et al.l 1)20101 ) ■ who 
used high-quality X-ray observations of Chandra and XMM- 
Newton to obtain temperature and density profiles of the hot 
gas. We find total mass distributions with an average mass 
of ~ 1.5 x 10 12 Mq and average circular velocity of ~ 515 
km/s at a radius of 25 kpc. 

The total circular velocity curves of our sample are all 
rising in the outer parts as a result of both an increasing 
temperature profile and a logarithmic pressure gradient that 
generally increases in magnitude outwards. Therefore the 
mass distributions of our sample are not isothermal, sug- 
gesting that we are probing the mass distribution sufficiently 
far out to observe the effects of the massive group/cluster 
haloes in which these galaxies reside. Increasing circular ve- 
locity curves appear to contradict work done in the past but 
are most likely a consequence of: i) Samples in the literature 
not restricted to only massive elliptical galaxies and/or ii) 
mass profiles not being probed sufficiently far out. 

The stellar components of our sample show a remark- 
able coherence in their surface-brightness profiles when 
scaled by the radius where d/i/dlog R = 2. We find however 
that the effective radius is a better proxy for the luminosity 
of the stellar component. The dark matter mass fraction is 
35-80% at 2R e and increases to about 80-90% at the out- 
ermost radii. 

Except in the case of NGC 1407, our determination of 
circular velocity curves from X-rays agrees well with previ- 
ous X-ray determinations. When comparing our circular ve- 
locity curves with those found by dynamical models in the 
literature, we find that the X-ray circular velocity curves are 
often lower, particularly in the central 10 kpc. This is prob- 
ably due to contributions from: i) Non-thermal pressures, 
ii) multiple-temperature components, iii) incomplete spatial 
coverage in the kinematic data used for the dynamical mod- 
els, and iv) insufficiently general mass profiles explored in 
the dynamical modelling. 

Looking at relations between the luminosity of the stel- 
lar component at R e , the total circular velocity at R e and 
the velocity dispersion of the environment, we find evidence 
for a dependence of the properties of central galaxies on their 
environment. 

The next steps in this project are to extend the analysis 
to a larger sample of X-ray bright elliptical galaxies to ex- 
amine the isothermality of their circular velocity curves, and 
to obtain more stringent constraints on how they compare 
with alternative methods of mass determinations, in par- 
ticular through the construction of more radially extended 
dynamical models. This will enable us to build up a clearer 
picture of the connection between properties of central X- 
ray bright elliptical galaxies and the environment in which 
they reside. 



X-ray bright elliptical galaxies 19 



ACKNOWLEDGEMENTS 

PD was supported by the DFG Cluster of Excellence "Origin 
and Structure of the Universe" . EC and IZ were supported 
by the DFG grant CH389/3-2 and OFN-17 program of the 
Russian Academy of Sciences. PD and IZ would like to thank 
the International Max Planck Research School (IMPRS) in 
Garching. 



REFERENCES 

Abadi M. G., Navarro J. F., Steinmetz M., 2006, MNRAS, 

365, 747 
Bender R., Dobereiner S., Mollenhoff C, 1988, A&AS, 74, 

385 
Brough S., Proctor R., Forbes D. A., Couch W. J., Collins 

C. A., Burke D. J., Mann R. G., 2007, MNRAS, 378, 1507 
Caon N., Capaccioli M., D'Onofrio M., 1994, A&AS, 106, 

199 
Cappellari M., Bacon R., Bureau M., Damen M. C, Davies 

R. L., de Zeeuw P. T., Emsellem E., Falcon- Barroso J., 

Krajnovic D., Kuntschner H., McDermid R. M., Peletier 

R. F., Sarzi M., van den Bosch R. C. E., van de Ven G., 

2006, MNRAS, 366, 1126 

Churazov E., Forman W., Jones C, Bohringer H., 2003, 
ApJ, 590, 225 

Churazov E., Forman W., Vikhlinin A., Tremaine S., Ger- 
hard O., Jones C, 2008, MNRAS, 388, 1062 

Churazov E., Tremaine S., Forman W., Gerhard O., Das 
P., Vikhlinin A., Jones C, Bohringer H., Gebhardt K., 
2010, MNRAS, 404, 1165 

Cowie L. L., Henriksen M., Mushotzky R., 1987, ApJ, 317, 
593 

de Lorenzi F., Gerhard O., Coccato L., Arnaboldi M., Ca- 
paccioli M., Douglas N. G., Freeman K. C, Kuijken K., 
Merrifield M. R., Napolitano N. R., Noordermeer E., Ro- 
manowsky A. J., Debattista V. P., 2009, MNRAS, 395, 
76 

de Lorenzi F., Gerhard O., Saglia R. P., Sambhus N., De- 
battista V. P., Pannella M., Mendez R. H., 2008, MNRAS, 
385, 1729 

Dejonghe H., de Bruyne V., Vauterin P., Zeilinger W. W., 
1996, A&A, 306, 363 

Drinkwater M. J., Gregg M. D., Colless M., 2001, ApJ, 548, 
L139 

Faber S. M., Wegner G., Burstein D., Davies R. L., Dressier 
A., Lynden-Bell D., Terlevich R. J., 1989, ApJS, 69, 763 

Forman W., Jones C, Churazov E., Markevitch M., Nulsen 
P., Vikhlinin A., Begelman M., Bohringer H., Eilek J., 
Heinz S., Kraft R., Owen F., Pahre M., 2007, ApJ, 665, 
1057 

Fukazawa Y., Botoya-Nonesa J. G., Pu J., Ohto A., 
Kawano N., 2006, ApJ, 636, 698 

Gavazzi G., Boselli A., Scodeggio M., Pierini D., Belsole 
E., 1999, MNRAS, 304, 595 

Gavazzi R., Treu T., Rhodes J. D., Koopmans L. V. E., 
Bolton A. S., Buries S., Massey R. J., Moustakas L. A., 

2007, ApJ, 667, 176 

Gebhardt K., Richstone D., Tremaine S., Lauer T. R., Ben- 
der R., Bower G., Dressier A., Faber S. M., Filippenko 
A. V., Green R., Grillmair C, Ho L. C, Kormendy J., 
Magorrian J., Pinkney J., 2003, ApJ, 583, 92 



Gebhardt K., Thomas J., 2009, ApJ, 700, 1690 
Gerhard O., Jeske G., Saglia R. P., Bender R., 1998, MN- 
RAS, 295, 197 
Gerhard O., Kronawitter A., Saglia R. P., Bender R., 2001, 

AJ, 121, 1936 
Gerhard O. E., 1994, in G. Contopoulos, N. K. Spyrou, 

& L. Vlahos ed., Galactic Dynamics and N-Body Simula- 
tions Vol. 433 of Lecture Notes in Physics, Berlin Springer 

Verlag, Elliptical galaxies, pp 191-274 
Hopkins P. F., Cox T. J., Dutta S. N., Hernquist L., Kor- 
mendy J., Lauer T. R., 2009, ApJS, 181, 135 
Humphrey P. J., Buote D. A., 2010, MNRAS, 403, 2143 
Humphrey P. J., Buote D. A., Brighenti F., Gebhardt K., 

Mathews W. G., 2008, ApJ, 683, 161 
Humphrey P. J., Buote D. A., Brighenti F., Gebhardt K., 

Mathews W. G., 2009, ApJ, 703, 1257 
Humphrey P. J., Buote D. A., Gastaldello F., Zappacosta 

L., Bullock J. S., Brighenti F., Mathews W. G., 2006, ApJ, 

646, 899 
Koopmans L. V. E., Bolton A., Treu T., Czoske O., Auger 

M. W., Barnabe M., Vegetti S., Gavazzi R., Moustakas 

L. A., Buries S., 2009, ApJ, 703, L51 
Kormendy J., 1977, ApJ, 218, 333 
Kormendy J., Fisher D. B., Cornell M. E., Bender R., 2009, 

ApJS, 182, 216 
Kronawitter A., Saglia R. P., Gerhard O., Bender R., 2000, 

A&AS, 144, 53 
Lauer T. R., Gebhardt K., Faber S. M., Richstone D., 

Tremaine S., Kormendy J., Aller M. C, Bender R., 

Dressier A., Filippenko A. V., Green R., Ho L. C, 2007, 

ApJ, 664, 226 
Liddle A. R., 2009, Annual Review of Nuclear and Particle 

Science, 59, 95 
Magorrian J., 1999, MNRAS, 302, 530 
Mahdavi A., Trcntham N., Tully R. B., 2005, AJ, 130, 1502 
Mandelbaum R., Seljak U., Cool R. J., Blanton M., Hirata 

C. M., Brinkmann J., 2006, MNRAS, 372, 758 
Matsushita K., Belsole E., Finoguenov A., Bohringer H., 

2002, A&A, 386, 77 
Merritt D., Tremblay B., 1994, AJ, 108, 514 
Naab T., Johansson P. H., Ostriker J. P., 2009, ApJ, 699, 

L178 
Nagino R., Matsushita K., 2009, A&A, 501, 157 
Nulsen P. E. J., Bohringer H., 1995, MNRAS, 274, 1093 
O'Sullivan E., Forbes D. A., Ponman T. J., 2001, MNRAS, 

328, 461 
Rix H., de Zeeuw P. T., Cretton N., van der Marel R. P., 

Carollo C. M., 1997, ApJ, 488, 702 
Romanowsky A. J., Strader J., Spitler L. R., Johnson R., 

Brodic J. P., Forbes D. A., Ponman T., 2009, AJ, 137, 

4956 
Saglia R. P., Kronawitter A., Gerhard O., Bender R., 2000, 

AJ, 119, 153 
Shen J., Gebhardt K., 2010, ApJ, 711, 484 
Spolaor M., Forbes D. A., Hau G. K. T., Proctor R. N., 

Brough S., 2008, MNRAS, 385, 667 
Thomas J., Saglia R. P., Bender R., Thomas D., Gebhardt 

K., Magorrian J., Corsini E. M., Wegner G., 2007, MN- 
RAS, 382, 657 
Thomas J., Saglia R. P., Bender R., Thomas D., Gebhardt 

K., Magorrian J., Corsini E. M., Wegner G., 2009, ApJ, 

691, 770 



20 P. Das et al. 

Thomas J., Saglia R. P., Bender R., Thomas D., Gebhardt 

K., Magorrian J., Richstone D., 2004, MNRAS, 353, 391 
Tonry J. L., Dressier A., Blakeslee J. P., Ajhar E. A., 

Fletcher A. B., Luppino G. A., Metzger M. R., Moore 

C. B., 2001, ApJ, 546, 681 
Trager S. C, Faber S. M., Worthey G., Gonzalez J. J., 

2000, AJ, 120, 165 
Tremblay B., Merritt D., 1996, AJ, 111, 2243 
Trentham N., Tully R. B., Mahdavi A., 2006, MNRAS, 369, 

1375 
Treu T., Koopmans L. V. E., 2004, ApJ, 611, 739 
van Dokkum P. G., Franx M., Kriek M., Holden B., Illing- 

worth G. D., Magee D., Bouwens R., Marchesini D., 

Quadri R., Rudnick G., Taylor E. N., Toft S., 2008, ApJ, 

677, L5 
Werner N., Zhuravleva I., Churazov E., Simioncscu A., 

Allen S. W., Forman W., Jones G, Kaastra J. S., 2009, 

MNRAS, 398, 23 
Zhang Z., Xu H., Wang Y., An T., Xu Y., Wu X., 2007, 

ApJ, 656, 805 



